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Abstract 

Motivated by the anomaly in the KARMEN experiment, we study new pos- 
\ sibilities for brane-world models of neutrino masses. We show that the KAR- 

MEN anomaly can be understood in the context of a six-dimensional brane- 
world model. The fine-tuning problem associated with the single-particle in- 
terpretation of the anomaly is thereby alleviated. Our model shows some 
' interesting properties of general interest for brane-world neutrino physics. In 

particular, the size of the brane-bulk mixing decouples from the scales relevant 
'"qJ f° r atmospheric and solar neutrino oscillations. 

43 ■ 

1 Introduction 

•rH , 

^ ' Two recent developments, namely the discovery of brane-world models [jl], ^, [ 

and the flexibility in choosing fundamental scales and sizes of additional dimensions 
||, |6|, have led to considerable activity exploring new options for particle phenomenology. 
Both developments have a natural place in string- or M-theory which leads to interesting 
connections between the leading candidate for a fundamental theory and particle physics. 

Much of the interesting phenomenology of brane-world models is associated with the 
Kaluza-Klein modes that originate from large, gravity-only additional dimensions. Par- 
ticularly the graviton and its Kaluza-Klein excitations are of interest in this context ||. 
Neutrino physics is another branch that is likely to be effected in brane-world models. 
This is because higher-dimensional bulk fermions lead to Kaluza-Klein towers of stan- 
dard model singlets that may be interpreted as sterile neutrinos |9|, [To], 11, 12]. Most of 




the work in brane-world neutrino physics to date has been focusing on rather specific 
models. Concretely, most of the explicit studies have been carried out for simple five- 
dimensional models with a special form of the action that, for example, does not include 
bulk mass terms (see, however, ref. ||, where Majorana mass terms have been considered 
in the context of orbifold models). 
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Brane-world models with large additional dimensions have led to a number of exciting 
predictions, such as the modification of gravity at small distances ||, that may be testable 
in the future. In this paper, we would like to pursue the idea that a signal originating from 
brane-world physics may have already been observed. Specifically, we will show that the 



anomaly in the KARMEN neutrino experiment [13, 14 1 can be explained in the context 
of brane-world models. This will also lead us to explore new possibilities for brane-world 
models of neutrino masses. That is, we will examine a number of new model building 
aspects in brane-world theories of neutrino masses. 

The KARMEN neutrino experiment studies neutral current interactions of neutrinos 
that originate from the decay chain ir + — ► /i + at rest. The KARMEN experiment has 
a unique time structure that allows one to study the number of neutrino events in the 
detector as a function of time after the production of the initial tt + . It is this distribution 
that shows an anomalous peak at a specific time after ir + production. The KARMEN 
anomaly can be explained in terms of a hypothetical X particle |l3| produced in a rare 
decay tt + — ► fi + X. It has also been shown [15] that the X particle can be interpreted 



as a sterile neutrino. While postulating such a particle leads to a satisfactory fit of the 
experimental data it has a theoretically unappealing feature. To match the observed peak, 
the X particle has to have a very specific mass mx that is fine-tuned to the mass difference 
m n — nia with an accuracy of approximately 10 -4 . That is, mx / (m w — m^) — 1 — 1.8 10~ 4 . 
There is no apparent theoretical explanation for such a coincidence of masses. 

Here is where brane-world physics might come into play. Suppose, that a brane- 
world model leads to a Kaluza-Klein tower of sterile neutrinos. Clearly, if the spacing is 
sufficiently small, it is relatively more probable for one of the particles in this tower to fall 
into the narrow mass range relevant for the KARMEN experiment than it is for a single 
particle. In the extreme case, where the spacing of the tower is of the order of the mass 
gap available for the X particle, the fine-tuning problem is solved completely. 

In this paper, we will show that a brane-world model explaining the KARMEN 
anomaly does indeed exist. An analysis of the basic requirements singles out two specific 
choices of dimensions and scales. The first one corresponds to a five-dimensional model 
with intermediate string scale and the scale of the fifth dimension in the keV range. The 
second one corresponds to a six-dimensional model with a TeV string scale and the scale 
of the two additional dimensions in the MeV range. We will focus on the second case 
and construct an explicit six-dimensional model that meets all the requirements. A major 
phenomenological constraint on this model is that the other Kaluza-Klein sterile neutrinos 
do not distort the KARMEN spectrum in a way inconsistent with the experiment. We 
will show that this can be avoided by making the spacing of the Kaluza-Klein tower suf- 
ficiently large. Unfortunately then, the fine-tuning problem cannot be solved completely 
in the way indicated above but is still improved by two orders of magnitude with respect 
to the single particle interpretation. At the same time, the effect of the Kaluza-Klein 
modes lighter than the X particle, provide a way of how our proposal could be experi- 
mentally distinguished from the single particle proposal. These particles contribute to the 
KARMEN spectrum at short times and, depending on the parameters of the model, the 
corresponding signal may be detectable in future experiments. 

Our six-dimensional model also shows a number of interesting properties that are of 
general importance for brane-world models of neutrino masses. In the simplest version of 
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the model with conserved lepton number, for example, we find three massless eigenstates 
that are generally nontrivial linear combinations of the electroweak eigenstates and the 
bulk neutrinos. The admixture of bulk neutrinos is controlled by a six-dimensional Dirac 
mass and can be made small or large while keeping the states massless. In contrast, in 
the models considered so far, massless modes were either not present or decoupled from 
the ordinary neutrinos. Furthermore, in those models, light neutrino masses and their 
mixing with bulk neutrinos are usually controlled by the same parameters. Our model 
shows that it is actually possible to decouple these two quantities. More generally, it 
is illustrated that the inclusion of all Lorenz-invariant (or even Lorenz non-invariant) 
mass terms in the higher-dimensional theory opens up new possibilities for brane-world 
neutrino phenomenology. 



2 The KARMEN anomaly and its brane— world interpretation 

The KARMEN experiment at the Rutherford Appleton Laboratory studies charged and 
neutral current interactions of neutrinos from the rr + — > /x + decay chain at rest. The 
unique feature of the neutrino flux is its time structure. The primary pions are produced 
in 0.5 /xs long pulses^ with a frequency of 50 Hz by the ISIS synchrotron proton beam. 
Due to the short pion lifetime, T n ~ 26 ns, the muon neutrinos from the ir + — ► /x + z^ decay 
represent the prompt component of the neutrino beam, which exhausts itself within a few 
tv after the pion pulse. Hence, in the time window (0.6-10.6) /xs after proton beam on 
target the neutrino beam is exclusively composed out of v e and Un that originate from 
the slow decay (r„ ~ 2.2 /xs) of the muons produced in the first 0.5 /xs. Consequently, the 
time spectrum of the number n of detector events in this time window is expected to be 
described by 

Here A is a time-independent^] rate constant depending, among other things, on the muon 
production rate, on the target cross section and on the detector efficiency. The constant 
B represents the event rate associated with a time-independent background. 

The KARMEN time-anomaly jL3L 14] manifests itself in a discrepancy between the 



expected time spectrum @j and the measured one. This discrepancy is mainly due to a 
peak at approximately 3.6 /xs after beam on target. Such a peak is clearly visible in Fig. |l]a. 
In this plot the time window (0.6- 10.6) /xs has been divided into 20 bins and the number 
of events for each bin (falling into the energy range from 10 to 36 MeV) has been indicated 
by a dot. In order to estimate the significance of the discrepancy and for later comparison 
with results based on our model we have performed a simple analysis of the data shown in 
Fig. [j]a. Specifically, we have carried out a x 2 fit of the data set in Fig. [l] assuming an event 
rate as in eq. dU) . Background measurements on a wider time interval fix the parameter B 
as indicated by the dashed line in Fig. |l|. The best fit results in Xmm = ^9 for 19 degrees 
of freedom, corresponding approximately to a 2a discrepancy. However, the significance 



1 More precisely, the 0.5 /is pulses are made of two shorter 100 ns long pulses separated in time by 330 ns 

2 The convolution of the muon exponential dec 
again an exponential with the same time constant 



2 The convolution of the muon exponential decay rate with the time spectrum of muon production is 
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Figure 1: Time distribution of events in the time window (0.6 - 10.6) fj,s after beam on 
target and in the energy interval (10-36)MeV. The time interval is divided in 20 bins. 
The dots represent the measured events. The dashed line indicates the measured back- 
ground. In Figure (a) the solid histogram shows the best fit of all bins except the 6th and 
the 7th (shadowed bins) using the event rate in eq. Q). In Figure (b) the signal associated 
with a Kaluza-Klein tower of states corresponding to X/R ~ 15MeV is added (see text). 



of the discrepancy is higher if one considers only the two bins around t = 3.6 fis. This 
is appropriate because an explanation of the anomaly around t = 3.6 //s by statistical 
fluctuation is strongly disfavored. The peak around t = 3.6 [is, already clearly present in 
the 1995 data |13]| , has, in fact, been confirmed by the subsequent KARMEN2 data [14]. 
Once the two bins around t = 3.6 /xs have been excluded, the remaining data is very well 
described by the distribution in eq. (jlj). The best fit leads to xLn = 14 with 17 de 

grees of 

freedom, corresponding to the solid line in Fig. |]a. According to such a fit, a total of 557 
events are expected in the two unfitted bins. As shown in [14], the measured 659 events 
represent a deviation at the level of approximately 4o\ Since the time-anomaly appears 
to be statistically significant and any attempt to explain its peak structure by systematic 
effects failed so far, it is worth to speculate about its physical origin. 

A simple description of the time spectrum can be achieved by assuming that the peak 
around t = 3.6 /us is associated with the decay of a slowly moving massive neutral particle 
"A" in the detector. Such a particle could, for example, be produced in the rare decay 
7r + —* fi + X pi]. The position of the peak then measures the time of flight of the particle 



and hence its velocity. According to the recent KARMEN fits [14| the velocity is given 
by fix — 0.0162. This velocity, in turn, determines the time the particle spends in the 
detector and, together with the time structure of the proton beam, the width of the signal. 
This width is found to be in nice agreement with the observed width of the peak. More 
generally, the slow particle interpretation is supported by a detailed analysis of the whole 
data set. Due to the very good time and spatial resolution of the detector, At < 2 ns, 
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Ax, Ay, Az < 15 cm, an analysis of the correlation between time and position of the 
events is possible. Such an analysis is in agreement with the slow particle hypothesis and 
provides the most precise measurement of the X particle velocity leading to the value 
given above 14]. 

Once the signal is attributed to the decay of a slow particle, a rare tt + — > [i + X decay 
represents the simplest possibility for the production of the particle. Since the time width 
of the anomaly is smaller than the muon lifetime, it is natural to assume that its production 
is not associated to muon decay (see however [|ll|). Moreover, the production in the 7r + 
decay allows to account for the slow, fixed velocity in terms of a small available phase 
space (unlike, for example, the production by the proton interactions with the target). 
Finally, the energy deposited by the particle in the detector is found to be less than about 
35 MeV [|14]] . This is what is expected from the decay of a particle produced in the decay 
7r + —* fJ- + X, whose mass mx has to be less than m n — m^ ~ 33.9 MeV. On the other hand, 
if the production process were tt + — > e + X, the slow X particle would have to have a mass 
mx — 137 MeV and would release more energy than measured (unless having a muon or 
another massive particle in the final state is for some reason favored). As for the decay, the 
simplest possibilities for the visible decay of a particle produced in the n + — > fi + X channel 
are X — > e + e~v and X — > jv, where v generically represents a neutral light particle. A 
visible decay mainly through X — > thereby producing mono-energetic photons is highly 
disfavored 114]. Consequently, in this paper we will consider the possibility that the X 
particle is produced in the decay ir + —* /i + X decay and the anomaly is mainly due to its 
visible decay X — > e + e~v. If the X particle lifetime is larger than 3.6 fis, as it will be in 
our case, the KARMEN data requires that 



BR(vr+ -> fi + X) • T vis = (1-2) HT 11 s~\ (2) 

where BR(-7r + — > /i + V) is the branching ratio of the decay ir + — > fi + X and r v j s is the 
partial width of the visible X decays. 



The first proposal for the nature of the hypothetical X particle was made in ]1l 
There, it was identified with a sterile neutrino that mixes with the Standard Model (SM) 
neutrinos. The same possibility was considered in [ll?], Supersymmetric scenarios were 



studied in [19, 20]. From a theoretical point of view, light sterile neutrinos are welcome if 
their lightness can be accounted for. Whereas the mass of the SM neutrinos is protected 
by the electroweak symmetry (and the absence of fundamental isovector Higgses), the 
SM symmetries alone cannot explain why sterile neutrinos should be lighter than, say, 
the Planck mass. It is therefore useful to consider sterile neutrinos in the context of a 
specific framework. In this paper, we will be dealing with sterile neutrinos that arise 
in brane-world models with large additional dimensions. Concretely, we consider the 
possibility that the KARMEN anomaly is due to a sterile neutrino that is part of a tower 
of Kaluza-Klein excitations associated with a SM-singlet fermion propagating in 4 + d 
dimensions. A possible explanation for small sterile neutrino masses is available if these 
models are viewed in the context of string theory. Then, bulk fermions arise naturally 
as superpartners of moduli fields. Perturbatively, moduli are flat directions of the theory, 
and, hence, their fermionic partners are massless before supersymmetry breaking. After 
including supersymmetry breaking effects those fermions will receive a mass related to 
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the supersymmetry breaking scale. Depending on the details, this might well provide a 
reason for small bulk fermion masses. The next sections will be devoted to a detailed 
study of such models, particularly in six dimensions. In the remainder of this section, we 
will describe some general features of our six-dimensional model from a phenomenological 
viewpoint and summarize its implications for the KARMEN experiment. 

Let us consider six-dimensional Dirac fermions as the origin of the sterile neutrinos. For 
simplicity we will consider only one such Dirac fermion here, although our later treatment 
will be more general. From a four-dimensional perspective, this field can be described 
by a Kaluza-Klein tower of sterile Dirac neutrinos v n labeled by two integer numbers 
n = (ni, 712). Assuming, for simplicity, that the two extra dimensions have equal size 2irR, 
the mass spectrum of the Kaluza-Klein tower is given by by 
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where n 2 = n\ + n 2 and the mass parameter \x originates from a six-dimensional Dirac 
mass term. Of course, to have a candidate for the X particle among these Kaluza-Klein 
fermions we have to require that the mass parameter \x is smaller than mass mx of the X 
particle. This requirement precisely corresponds to the problem of keeping sterile neutrinos 
light. With the above string theory explanation for this in mind, we assume that \i < mx- 
Then, we may have a specific sterile neutrino X = u nx in the Kaluza-Klein tower that we 
would like to associate with the KARMEN anomaly. From eq. (||) we can compute the 
density of states p(mx) at the scale mx, defined as the average number of states dn per 
mass interval dm. Neglecting the mass \i for an order of magnitude estimate, this density 
can be computed from eq. @ to be | 

p(mx) ^ 2irR 2 mx ■ (4) 

In other words, the average separation between two subsequent states at mx in a mass- 
ordered list is given by 

KK p(mx) 2nR 2 mx ^ ^ 

The last equation highlights an interesting feature of our framework. Besides being able to 
account for the existence of light sterile neutrinos, it allows us to alleviate the fine-tuning 
problem associated with the slow particle interpretation of the KARMEN anomaly. The 
fine-tuning problem arises because the slowness of the X particle requires its mass to be 
only slightly lower than the difference m n — m^ ~ 33.9 MeV. More precisely, in order to 
be produced with a velocity fix — 0.0162, the X particle mass mx must be 

m x = m n - m^ - 5m x , (7) 

where 

5m x ^ 0x ("V ~ "I/Jtt^- - 6keV. (8) 
2m M 



3 In our simplified model the n-th state is eight times degenerate if m ^ 0, ni =^ 0, n\ 7^ 712- In a 
general situation where the size of the two extra dimensions are different and all possible six-dimensional 
mass terms are included, the degeneracy is lifted. 
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This means that mx is fine tuned to the — mass difference with a precision of 
1.8 x 10 -4 ! On the other hand the fine tuning problem is significantly alleviated in the 
approach we propose, since the X particle is just one of the many excitation of the six- 
dimensional fermion. It is therefore relatively more probable that one out of the many 
states has a mass that falls into the critical range. 

In principle, this framework allows us to completely eliminate the fine tuning problem. 
In fact, the problem disappears if the average separation between subsequent states is 
about twice the mass range available for the X particle, that is, 5mx = m n — — mx ~ 
Awkk/2- Prom eq. (|5|) this happens if 

— = (2Trm x Am KK ) 1/2 ~ (4irm x 5m x ) 1/2 ^ 1.6 MeV. (9) 
R 

In this case, the production of the particles heavier than X in the ir + — > [i + v n decay would 
be kinematically forbidden^. However, the states lighter than mx, being faster than the 
X particle, modify the time spectrum in the window (0.6-3.6) fj,s. In order compute the 
effect of those lighter states on the time spectrum, one has to take into account that 
they are produced more copiously due to the larger phase space available. However, their 
signal is shorter because they spend less time inside the detector. A simulation taking into 
account the detailed time-shape of the signal associated with the lighter particles shows 
that the KARMEN data is in very good agreement with the mass spectrum associated with 
an inverse radius of about 15 MeV, namely about 10 time larger than the natural value 
in @0. Fig. |l|b shows the corresponding best fit (solid line). Once the mass spectrum 
has been fixed by the choice of 1/R, the overall normalization of the sterile neutrinos 
signal and the constant A in eq. (jl]) are fitted to the data. The best fit shown in Figure 
has a Xmm = 15-8 for 18 degrees of freedom. The fine tuning problem is alleviated. The 
maximum density of states at the mx scale is about 40 times larger than (m n — m^) _1 , 
making this framework about 40 times more natural than the standard framework with a 
single particle. However, a significant fine tuning is still needed. 

The discussion above points out a peculiar feature of our proposal. The X particle 
comes with a number of particles lighter than mx whose number depends on the radius 
of the additional dimensions. For large values of the inverse radius, the lighter states are 
separated by a large mass gap. In this case our model is indistinguishable from models with 
a single sterile neutrino as far as the KARMEN experiment is concerned. However, from 
the theoretical point of view, it still provides a well motivated scenario for the existence 
of a sterile neutrino with a mass of 33.9 MeV which, at the same time, can account for the 
standard oscillation phenomenology (see Section Q). On the other hand, if the average mass 
gap Am-KK between the states is smaller than about 5 MeV, the particles lighter than X 
could give a small contribution to the early KARMEN time spectrum. For example, with 
the parameters used for the fit shown in Fig. |l]b, corresponding to an average separation 
of about 1 MeV between states, one gets a small contribution to the number of events in 
the first time bin. Therefore, a detailed study of the time distribution at early times after 
the end of the proton pulse could test the model if the mass separation Atti-kk happens 

4 This assumes that there is no local inhomogeneity in the density of states, so that all states heavier 
than X are indeed heavier by a margin larger than the mass gap 8mx of the X particle. 

5 The simulation assumes that the degeneracy mentioned in footnote His completely lifted. 
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not to be too large. An upgraded detector with tracking capability devoted to the study 
of the anomaly could certainly explore this possibility. Needless to say, the detection of a 
signal at early times would represent a strong hint for the existence of a brane-world. 



3 Scales and couplings 

In this section, we would like to analyze the scales and couplings of the brane-world theory 
that are required for a solution of the KARMEN anomaly along the lines explained above. 
This analysis will provide us with an overview over the various models that may be suitable 
for our purpose. In the following section, we will pick one of the models found in this way 
and develop it in detail. 

We start with a ten-dimensional string theory with string scale M s . The six internal 
dimensions are schematically divided into two groups. The first group of d dimensions are 
the ones that give rise to Kaluza-Klein neutrinos potentially relevant for the KARMEN 
experiment. The other 6 — d dimensions are those without any direct relevance for the 
experiment. For example, a particular dimension is not relevant in this sense if its inverse 
size is much larger than the mass of the X particle. For simplicity we assume that the first 
d dimensions have the same size characterized by the radius R. This is likewise assumed 
for the remaining 6 — d dimensions where the corresponding radius is called p. Given 
this setup, there are three basic requirements to be satisfied. First, the four-dimensional 
Planck constant Mp\ specified by 

(ir) = 16 ^ RM s) d (^pMs) G - d (10) 

should have the correct value. Except for this general condition we have two more con- 
straints that are specific to our problem. The first one is related to the probability of 
having a Kaluza-Klein particle with mass in the relevant range for the slow-particle inter- 
pretation of the KARMEN experiment. We would like to significantly reduce the amount 
of fine tuning required. This means that the parameter A defined as 

— = p(m x )Ara K K (11) 

should be significantly smaller than O (l0 4 ) , corresponding to the fine-tuning of the single- 
particle solution. We remind the reader that p(mx) is the density of Kaluza-Klein states 
at mx = 33.9 MeV and Attikk — 6keV is the mass gap between the X particle mass and 
the n-p mass difference. Note that, from the distinction of dimensions made above, only 
the Kaluza-Klein modes associated to d dimensions are relevant for this density. For iV 
Dirac spinors in 4 + d dimensions with a Dirac mass p it is given by 

jo[d/2] d/2 

P(mx) = T{d/2 + l) NRd ^™l ~ ^) d/2 - 1 ■ (12) 

The final requirement originates from the constraint @ on the branching ratio and the 
visible decay width of the X particle. This constraint implies that there should be a non- 
negligible mixing between the X particle and the left-handed neutrinos. In the models 
under consideration such a mixing is generated by a mass term between brane and bulk 



S 





(iVAy 2 )- 1 /^ 
M s 


(iVA)- 1 /^ 
iT 1 


d-S -IS 

(NA) d{f>-d)yd{6-d) X 


d= 1 


4.5 x 10 Y GeV 


12keV 


7.9 x 10 B GeV 


d = 2 


1.2 TeV 


1.6 MeV 


461 MeV 


d = 3 


23 GeV 


5.6 MeV 


28keV 



Table 1: Shown are scales and dimensions relevant for an explanation of the KARMEN 
experiment. We present typical values of the string scale M S; the inverse radius R~ l of the 
d additional dimensions associated with relevant Kaluza-Klein neutrinos and the inverse 
radius p^ 1 of the remaining 6 — d dimensional internal space. 

fields that resides on the brane. On dimensional grounds the associated mass parameter 
m can be written as 

hv = hmt 
m {2nRM s ) d / 2 h t {2itRM s ) d l 2 ' 1 ' 

where h is a brane-bulk Yukawa coupling, v is the the standard model Higgs VEV and 
ht and mt are the top Yukawa-coupling and mass. Then, the mixing matrix element U v x 
between the left-handed neutrinos and the X particle specified by 

U vX = — (14) 

mx 

should not be too small. Typically, from eq. @j one has \U u x\ — 10~ 3 . 

Altogether, we now have three constraints for our three fundamental parameters M s , 
R and p. All other quantities appearing in those constraint are either fixed or in a well- 
defined range. We can therefore simply solve for M s , R and p as a function of the number 
of dimensions d. The fact that we can basically determine all fundamental scales is quite 
remarkable and shows how constrained the problem is. A priori, it is not clear at all that a 
solution with sensible values of the scales exists. However, as Table [l] shows, this is indeed 
the case. The quantity y appearing in the table is defined as 

y= Tm^h t - (15) 

Let us discuss the various cases shown in table |l[ For one additional dimensions d = 1, 
we have a model with intermediate string scale, one large additional dimension in the keV 
range and the size of the other five dimensions close to the string scale. For two additional 
dimensions, d = 2, the situation is quite different. We have a model with the string scale 
in the TeV range, two large additional dimensions in the MeV range and the scale of the 
other four internal dimensions being about two to three orders of magnitude larger. In 
the case of three additional dimensions, d = 3, the typical string scale is already quite low 
and one needs a large factor NAy 2 to elevate it above the required TeV limit. Therefore, 
this case only represents a marginal possibility^. Finally, one can show that models with 

6 Moreover, unlike for the other two cases, the scale p _1 is smaller than the mass mx of the X particle. 
Therefore, if the Kaluza-Klein neutrinos propagated in those dimensions they will be relevant for the 
experiment, contrary to our initial definition of those dimensions. Therefore, one has to assume that this 
is not the case. Such an assumption is probably not very appealing from the perspective of model-building. 
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more than three additional dimensions, d > 3, lead to a string scale below TeV and are, 
therefore, not viable in our context. 

In summary, we have isolated two cases which may lead to an explanation of the KAR- 
MEN anomaly in the context of a brane-world model. The first case represents an effective 
five-dimensional model with intermediate string scale and the remaining five dimensions 
being close to the string scale. The second one is a model with a TeV string scale, two 
large dimensions in the MeV -1 range and the scale of the remaining four dimensions 
being two to three orders of magnitude larger. Therefore, this model is effectively six- 
dimensional in an intermediate energy range. We stress that the existence of these two 
possibilities, particularly the one with a TeV string scale, was by no means guaranteed 
and is due to a rather fortunate conspiracy of the various mass scales in the problem. 
While explicit examples can be constructed in both cases, in this paper we will focus on 
the six-dimensional case with a TeV string scale. 



4 An explicit six dimensional example 

In this section, we would like to present an explicit brane-world model that realizes the 
ideas outlined above. 

It is clear from the above discussion of scales that it is sufficient for our purpose to con- 
struct an effective six-dimensional brane-world model, valid for energies below p~ l . The 
Kaluza-Klein modes associated to the remaining four internal dimensions are too heavy to 
be relevant within our context. More specifically, we would like to couple a six-dimensional 
bulk with right-handed neutrinos and a three-brane that carries the standard-model fields. 
Such a model should then be analyzed in detail with respect to its consequences for the 
KARMEN anomaly. Furthermore, as we will see, the model illustrates some general fea- 
tures that arise in brane-world models of neutrino masses which have not been considered 
so far. 

The six-dimensional bulk action of the model is specified by 

Sbuik = J d 4 xd 2 y l^^idA^i - (hij*li*rj + h.c.)] . (16) 

Focusing on the relevant Yukawa couplings between bulk and brane fields, the brane action 
reads 



'S'branc — / d X 

J{y=0} 



^f a7 L^ + h.c. 

M s 



(17) 



We use coordinates x A with indices A,B, - - = 0, • • • ,5 for the total six-dimensional 
space-time. Four-dimensional coordinates x M are indexed by - = 0,1,2,3. The 

remaining two coordinates are denoted by y a , where a,/3, ■■■ = 1,2. The signature of 
our metric is "mostly minus". We have introduced N six-dimensional Dirac fermions 
Vl/j, where I , J, ■ ■ ■ = 1, ■ ■ ■ , N, in the bulk. Their left- and right-handed components are 
defined in the usual way by ^r/l,i = (lilV) 1 ^/ where Ta are the six-dimensional gamma 
matrices and = Tq ■ ■ ■ T5. The 4 + 2 decomposition of these gamma matrices can be 
written in the form Tj = {7^ (g> 12,275 <£> ff a }- Here 7^ are the four-dimensional Dirac 
matrices and 75 = — 270717273- The two-dimensional Dirac-matrices a a can be identified 
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with the first two Pauli matrices. Defining the two-dimensional chirality operator by 
03 = — iaia2 we have the relation Ty = 75 ® 03 between six-, four- and two-dimensional 
chiralities. For later purposes we also define six-dimensional charge conjugation by = 
C~ 1 *f r J\ where the charge conjugation matrix C is specified by the relations (Xa) T = 
— CT^C" 1 , as is usual in even dimensions. Furthermore, we have introduced Dirac mass 
terms with associated mass matrix jjlij in the bulk. By a suitable redefinition of the 
bulk fermions we can diagonalize this mass matrix. In the following, we will use this 
diagonalized form 

fiu = diag(//i, ■ ■ ■ ,/ijv) (18) 

with real, positive eigenvalues fij. We also choose the flavor basis for the lepton doublets 
Li such that the charged lepton Yukawa couplings are diagonal. In order to be able to 
couple the bulk spinors to brane field we decompose each six-dimensional field into two 
four-dimensional Dirac spinors *S> a i where a = +, — indicates the internal two-dimensional 
chirality of the component. The three-brane is taken to be located at y = in the internal 
space and carries, among the other standard model fields, the lepton doublets Li, where 
i,j,--- = e, fi, r, and the Higgs doublet H. These fields have a Yukawa coupling to the two 
components \l/ a j of the bulk spinors introduced above where the dimensionless coupling 
constants are denoted by h a n. Note that, from the dimensionality of the bulk spinors, the 
corresponding operator are suppressed by one power of the string scale M s . 

In summary, our model consists of N "right-handed neutrino" bulk spinors with a 
Dirac mass in six dimensions and Yukawa couplings to the lepton doublets located on the 
three-brane. Our model, as stands, has a lepton number U(l) symmetry with Li and *$>i 
each carrying one unit of charge. We would like to impose this U(l) symmetry (or an 
appropriate discrete subgroup thereof) on our model. This forbids the other mass terms 
and couplings that we could have written in our action. In particular, it forbids the bulk 
Majorana mass terms *$> c i*$>j and ^ c [T^j which would be allowed by six-dimensional 
Lorenz invariance. Furthermore, it forbids the brane-bulk coupling ^ c a iLiH that would 
be allowed from four-dimensional Lorenz invariance. However, we stress that, at this 
stage, nothing forbids the bulk Dirac mass term. It has, therefore, been included in the 
above action. 

Next, we would like to work out the four-dimensional effective action of our model. We 
compactify the additional dimensions on a two-dimensional torus that, for simplicity, is 
taken to be rectangular and with equal radii R in both direction. With the Higgs vacuum 
expectation value v, we introduce the mass parameters 

haliV . . 

maU = 2^RM S ■ (19) 

We expand the bulk spinors into Kaluza-Klein modes as 

2/) = ^E ex P ) (2°) 

where n = (ni,ri2). As we did with with the spinor ^/ before, we decompose each of its 
Kaluza-Klein modes ^>i n into two four-dimensional Dirac spinors ^ a in, where a = +,—. 
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Each of these four-dimensional Dirac spinors is decomposed further into Weyl spinors 
according to 



^aln 
Vain 



(21) 



The charge conjugated Weyl spinor £ c is defined as £ c = e£* where e is the two-dimensional 
epsilon-symbol. As a result, for each bulk spinor and each fixed mode number n we 
obtain two four-dimensional Dirac spinors ^f a j n or, equivalently, four left-handed Weyl 
spinors iaim Vain, where a = +, — . With these definitions, the mass terms in the four- 
dimensional effective Lagrangian read 

£bulk = ^2 i m nt+InV-In ~ m* n ^- In r] +In + fll(£+I n r] + l n + £-InV-In) + h.C.] (22) 
£branc = ^ [ m +I^+In + m_/i£_/ n ] + h.C. . (23) 

nez 2 

Here, the Kaluza-Klein masses m n are defined by 

mi + n 2 



R 



(24) 



To interpret the structure of masses given above it is useful to diagonalize the bulk La- 
grangian. Since for each mode number n and bulk flavor / we have two degenerate Dirac 
spinors, there is a degree of arbitrariness in the choice of the eigenstates. To keep the 
brane Lagrangian unchanged, we rotate the spinors r\ only. We therefore introduce the 
linear combinations 



/ ViV+in + m n ri_ In 

V +i n ~- 



i ) V —In I 

^H 2 j + \m n \ 2 \J^ 2 i + \m n \ 2 

which allows us to write the bulk Lagrangian in the form 



VlV-ln ~ rrf n r) +In 



(25) 



£bulk = ^2 y At/ + \m n \ 2 [(,+InV+In + Z-InV'-In + h - c -] 

n£l? 

From this we read off the following mass matrix 



i+n 
M = V'+n 

C— n 
V'-n 



( o 

rh + 




mi 












\J fJ 2 + \ m n\ 2 



yj '/J 2 + \m n \ 2 








- T 

m_ 









2,1 |2 
A 1 + \ m n\ 









2,1 |2 
A 1 + \ m n\ 



(26) 
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(27) 
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Can — (Cain ■ • • £i 



aNn) 



Van 



where for ease of notation we have used v = (v e ,Vn 
(Vain • • • v' a Nn) T ana - f 1 = diag(^i ... /i at). Given this matrix, it is easy to show that we 
have three exactly massless Weyl fermions VQi given by 



rn 



alj 



n^a \[$ + \m n \ 2 



Z V aln 



(28) 



Since in our context all the mass parameters are several orders of magnitude above the 
eV scale, it is natural to attribute the solar and atmospheric neutrino phenomenology 
to these massless eigenstates. The occurrence of massless eigenstates is due to the U(l) 
symmetry we have imposed. We will discuss below how a small breaking of this symmetry 
can generate the small masses necessary to account for the standard phenomenology of 
neutrino oscillations. As long as the states ^oi are massless, the matrix N in eq. ( p8|) 
can be determined only up to a unitary transformation among the massless eigenstates. 
However, N will be fixed (up to phase rotations) once small masses for the v^i have been 
generated. After integrating over the Kaluza-Klein modes up to the cut-off M s one finds 



in 



ali m alj 



a.l 



- 2 +R 2 



In 



M 2 



^ + 1/R 2 



(29) 



We are intereste d in a situa tion where the mass matrix has a hierarchical structure sat 
isfying \m a n\/ 'J /j, 2 + \m n \ 2 <^ 1 for all states. In this case (N^N)^ 
approximately given by a unitary matrix U, that is 



5ij and N is 



N ~ U 



(30) 



An exact diagonalization of the mass matrix A4 shows that the massive eigenstates can 
be organized in two Dirac spinors for each bulk flavor index I and each mode number 
n. Assuming a hierarchical structure of the mass matrix also simplifies the process of 



finding such eigenstates since it implies that the mass matrix (27) can be diagonalized 



perturbatively. While this reproduces eq. (28) for the massless states, it tells us that 
for each bulk flavor index / and each mode number n the two Dirac spinor eigenmodes 
iKin^ainf are given approximately by 



v aln ~ 
Vain ^ 



Cain 
Vain + 



fJ-j + \m n \ 2 



-Vi 



(31a) 
(31b) 



m. 



where a = +,— . This equation holds up to corrections of order \m a n\ 2 / \\x\ + 
Furthermore, in the same approximation, these two spinors have a degenerate Dirac mass 



M Ir 



+ \m n \ 2 



(32) 



Given this information we can invert the formula ( |28| ) and express the weak eigenstates vi 
in terms of mass eigenstates VQi, v a i n as 



Vi ~ UijV j + ^2 



m 



al i 



n,I,a \l ~\~ l^nP 



-Vain 



(33) 
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up to small terms of order \rh a ji\ 2 / (fi 2 + |m n | 2 ). Eq. (|33| ) is the starting point for studying 
the phenomenology of our model. Due to our convention for the flavor basis of the left- 
handed neutrinos (leading to diagonal charged lepton Yukawa couplings), the V{ are the 
neutrino "flavor eigenstates" . They are expressed as a superposition of three light Majo- 
rana mass eigenstates (y$i) and a series of left-handed components of heavy Dirac mass 
eigenstates (v a ln)- I n the limit \rh a n\ 2 / (/i 2 + |m ra | 2 ) <C 1, the "heavy" contribution to the 
states Vi is small, so that the flavor eigenstates are mainly light. Therefore, on one hand 
the model can potentially accommodate the standard oscillation phenomenology, which 
mainly involves the three light Majorana eigenstates z^oi an d the 3x3 mixing matrix U. 
On the other hand, as we will see, the small admixture of the Kaluza-Klein states account 
for the KARMEN anomaly. 

Let us discuss the results that we have obtained so far from the general perspective 
of brane-world models for neutrino masses. We have found three exactly massless Weyl 
fermions as given by eq. (28). In general they are superpositions of the electroweak eigen- 



states Vi on the brane and the Kaluza-Klein states v( aln . The relative weight of these two 
contributions is controlled by the Dirac masses fij. For vanishing Dirac masses /U/ — > the 
right hand side of eq. (^) is dominated by the terms corresponding to the Kaluza-Klein 
zero mode r]' aIQ . Hence, in this limit, the appearance of massless modes is not surprising. 
The massless states ( |28| ) are simply linear combinations of the bulk zero modes in this case. 
The situation becomes less trivial once we switch on fij. Then, there are no massless bulk 
modes any more since, from eq. (32), the Kaluza-Klein spectrum is bounded from below 



by To summarize, despite the introduction of a bulk Dirac mass term, the symmetries 
of the higher-dimensional theory (Lorenz invariance and lepton number) lead to massless 
eigenstates which are non-trivial linear combinations of the brane and the bulk fields. We 
believe this can be of considerable relevance for brane-world models of neutrino masses, 
in general. A final point concerns the order of magnitude of the bulk Dirac mass \x. As 
stands, the natural value of this mass is probably the string scale M s . Of course, in this 
case, the bulk neutrinos would be too heavy to be relevant for any neutrino physics. If 
the string scale is of order TeV some bulk states might receive smaller masses of order, 
say, MeV due to small couplings. In addition, if the bulk fermions are interpreted as 
superpartners of string moduli they receive masses only after supersymmetry breaking, 
which may also account for their smallness. 

So far, we have simply assumed that the massless states acquire a small mass through a 
breaking of the U(l) symmetry. However, we have not explicitly shown how light neutrino 
masses and mixings needed for the standard neutrino phenomenology can be generated. 
While this is certainly appropriate for our main purpose it is, of course, important to show 
how small but non-vanishing neutrino masses can be incorporated into the model. This 
is what we will do in the remainder of the section. 

In order to generate neutrino masses we need to break the lepton number symmetry 
that we have imposed on our model. Rather than presenting a complete model of how 
this can be realized spontaneously, we simply parameterize this breaking by adding the 
following Majorana mass terms to the bulk action (|l6|) 



-S'buik.L = / d 4 xd 2 y 



l -M\f^j - ^M^^T^j + h.c. 



(34) 
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(A) (S) 
with Mjj antisymmetric and Mjj symmetric. They violate lepton number by two units. 



Furthermore, we have to amend the brane action (|17j) by 

(35) 



' brane, L 



-^^ + h.c. 



d 4 x 

l{y=0} 

which violates lepton number by two units as well. Recall that ^> a i, where a = +, — are 
the two four-dimensional Dirac spinors contained in the six-dimensional bulk spinors ^j. 
As before, after electroweak symmetry breaking, we introduce the mass parameters 

laliV toci \ 

maR = 2^RM S ■ (36) 

Using the expansion (^) of the bulk fermions we then find the additional four-dimensional 
terms 

£buik,L = ^2 \M u r] +In r)-j- n - M}j€ +Jn €-j-n + h.c. (37) 

nSZ 2 

£brane,L = ^ [m+IiC+In + rh-Ii€-In] + h.C. (38) 

nez 2 

where the Major ana mass matrix M is defined by M = + . Now our total four- 
dimensional effective action consists of the original lepton number preserving parts (|22|), 
(|23|) and the above lepton number violating parts (|37|), (j38j). It would clearly be interesting 
to thoroughly investigate the neutrino phenomenology resulting from this action that 
includes all obvious bulk mass terms allowed by higher-dimensional Lorenz invariance. 
We hope to return to this problem in a future publication. For the purpose of the present 
paper, we would merely like to check the order of magnitude of masses that are induced. 
To do this, we focus on the case of one flavor in the bulk as well as on the brane. Assuming 
that m a /fi <C 1 and rh a j '/■* <C 1 we can apply the see-saw formula to find the light neutrino 
masses. After summing over all Kaluza-Klein states with cut-off M s we find for the mass 
mo of the previously massless states i/oi 

rriQ = {—rh-rh-H — rh-rh+M + rh—rh+M* — m+m-f^u) 

M 2 \ 
+ 2KR 2 \n— . . 9 S . _ . (39) 



,/x 2 + |M| 2 /x 2 + \M\ 2 + \/R 2 i 

As expected, this vanishes if we restore lepton number by setting fh a = M = 0. The term 
in the first parenthesis together with the first term in the second parenthesis correspond 
to what one normally expects for a see-saw suppressed neutrino mass. In fact, these terms 
describe the contribution from the lightest Kaluza-Klein state with mode number n = 0. 
However, in addition we have a large number of higher Kaluza-Klein states each of which 
contributes to the light neutrino mass. The net effect of all these states is encoded in the 
last term in eq. (|39j). This term is proportional to R 2 rather than the explicit mass scales 
li or M. Hence it corresponds to a suppression by the mass scale 1/R associated to the size 
of the additional dimensions. As we will see, in our context, this suppression is sufficient 
to obtain reasonable small neutrino masses under plausible assumptions. Eq. ( p9| ) also 
illustrates the above mentioned decoupling of neutrino masses and neutrino-bulk mixing. 
The former depend in particular on the lepton number violating quantities m± and M. 
As long as these quantities are sufficiently small the light neutrinos are still approximately 
given by eq. (|28|) and, hence, they depend on lepton number conserving parameters only. 
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5 Quantitative analysis of the model 



We would now like to demonstrate more quantitatively that the six-dimensional model 
presented in the previous section can explain the KARMEN anomaly and is, at the same 
time, compatible with other phenomenological constraints. 

Let us start by explaining the requirements on our model that follow from the KAR- 
MEN anomaly. This amounts to specifying the details of the fit that we have presented 
in Section Q and work out some of its consequences. Let us start with the constraint (|2|) 
on the branching ratio and on the decay width of the X particle that our model has to 
satisfy. We should express this constraint in terms of the parameters of our model. The 
X particle corresponds to a mode nx such that 



m X = ^//i 2 + |m nx | 2 . (40) 

As we saw in the previous Section, two almost degenerate Dirac mass eigenstates with 
opposite internal chirality are associated to this mode nx, namely (^+ nx , v +nx ) T and 
(vt n , v- nx ) T . The mixing U? x , where a = +, — between the SM flavor eigenstates i/j, 
i = e, fx, r, and these two states can be read off from eq. (^) as 

\m = . &d = (4i) 

V + K*! 2 mX 

However, the two states are not completely degenerate but split at higher order in the 
see-saw approximation. For some of the parameter space this splitting will be smaller 
than the X particle mass gap 5mx — 6keV. In this case, both states contribute to the 
KARMEN anomaly and can effectively be accounted for by a single state X with mass 
parameter 

m; = \m +i \ 2 + |m_i| 2 . (42) 

If, on the other hand, the splitting is larger than brax — 6keV one mass eigenstate will 
be above the threshold for production in the KARMEN experiment and the anomaly is 
due to the other state. We will call the corresponding mass parameter rh as well, but it 
is now a function of rh±i generally different from the one given above. For both cases, we 
define the mixing matrix elements by 

\U iX \ = ^ ■ (43) 
m x 

Through those mixings, the mass parameters m, control the X particle production and 
decay and therefore the number of anomalous events in the KARMEN experiment. The 
branching ratio for the decay ir + — > ll + X is given by 

\( 7TV^ Tup 1 Trip 1 1 

BR(tt + - u+X) = |CV| 2 \ f % ^ n f * 0.028\U, X \ 2 , (44) 

where 

A(a, b, c) = [a(b + c) - (6 - c) 2 ] V 'a 2 + b 2 + c 2 - lab - 2ac - 26c. (45) 
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Possible contributions from bulk dynamics to the visible decay width are highly suppressed 
compared to those from standard W and Z exchange. The impact of bulk dynamics on 
the total width can in principle be non-negligible but it would not change drastically 
the lifetime, which means that it would not affect our conclusions at all. Neglecting the 
sub-dominant X — > decay, the visible decay width is given by 



r ■ ~ 

L vis — 



F(X -» e + e~v) 

5 



m x 



m. 



\Uex\ 2 



-7 + Sw + Isyy I + (\U^ X \ 2 + \U tX \ 2 ) ( ~ 



_-l 



890sec _1 |?7 e x| 2 + 195 sec" 1 (| f/ M j^ | 2 + \U rX \ 2 ). (46) 



Putting these results together, the constraint (0) can be expressed in terms of the mixings 
\Uix\ as follows 



(1-2) 10 



-ii 



BR(tt+ -> ^+X)r vis sec 



2h\U llX \ 2 \Uex\ 2 + 5.5|C/ M x| 2 (|^x| 2 + \U T 



X 



(47) 



If we make the plausible assumption m 2 <C m 2 ,m 2 , it follows that the first term on the 



right hand side of the previous equation is sub-dominant. Eq. ([471) then represents a 
constraint on |£/^x| an d |t/"rx|> or > equivalently on rh^ and m r which reads 



[(40-47) keV] z 



(48) 



This constraint relates and rh T and allows us to compute BR(-7r + — * n + X) and r v ; s as 
a function of the ratio 



r = rh^/rh T 



(49) 



The result is shown in Fig. ^ for the range 10 -4 < m^/m T < 1 which corresponds to 
the plausible situation < rh T . The width of the band in these figures corresponds 
to the uncertainty on the right hand side of eq. (|48|) . If the "atmospheric" mixing angle 
measured by SuperKamiokande originated in the diagonalization of the charged lepton 
mass matrix, we would expect rh^/m T = O (1). This is because the brane-bulk Yukawa 
coupling was written in a basis where the charged lepton mass matrix is diagonal. However, 
w«/w T = O (1) gives a too high branching ratio BR(7r + — > H + X) according to the present 
PSI upper limit (solid horizontal line in Fig. |^a). The new expected limit is also shown 
(dashed horizontal line in Fig. ||a). The latter would imply rh^/rh T < 10~ 2 . Therefore, 
we focus on the case where m^/rh T <C 1. Furthermore, a lower constraint on the m M /m T 
range is given by the requirement that rh T is not too large. If rh T were of the order of 
the Dirac mass fi or more, the r neutrino would predominantly consist of a heavy mass 
eigenstate, which would be phenomenologically problematic. Consequently, r should be 
preferably in the range 



10 



- 4 to lO" 2 



(50) 
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Figure 2: Branching ratio BR(ir + — > (a,), and inverse visible width T v i s , (b), as 

functions of the ratio r = rh^/m T . 

Then, the mass parameters rhi of our model as well as the branching ratio and the decay 
width of the X particle can be expressed in terms of r as 

1/2 



m M = (4.0~4.7)keV(^ I ^ f J (51a) 

/10-2\V2 

m T = (0.4 - 0.47) MeV (51b) 



r 



BR(vr + -► ll + X) = (3.8-5.3) 10~ 10 [j^J (51c) 
r x <r v - is 1 = (27-38)sec (^) . (51d) 

It remains to fix the radius R of the two additional dimensions. Of course, 1/R should 
be smaller than the mass mx — 33.9 MeV of the X particle. Furthermore, as we have 
seen in Section ^, a value of 1/R ~ 1.6 MeV leads to a complete solution of the fine- 
tuning problem. However, as explained there, we need to choose a somewhat larger value 
1/R> 15 MeV to avoid significant contributions from modes lighter than the X particle 
to the KARMEN spectrum. In summary, we therefore require 

15 MeV < 1/R < 33.9 MeV . (52) 

The fine-tuning is smallest at the lower end of this range which is why we have chosen 
1/R ~ 15 MeV in the fit performed in Section ^. Finally, we have to fix the order of 
magnitude of our Dirac mass /x. To avoid too large admixture of the lightest Kaluza- 
Klein mode we have to require that \x 3> rh T . In addition, of course, \i has to be smaller 
than the X particle mass. Together, this leads to 

rh T < fj, < 33.9 MeV . (53) 

We remark that this allows for the plausible value fi = 0(1/R) for the Dirac mass. Finally, 
we should check the assumption made in the previous section about the unitarity of the 
matrix iV and the smallness of the heavy states contribution to the SM neutrinos. Both 
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issues are related to the square of the normalization matrix N^N defined in eq. (26) which 
should be close to unity. For our favorite values 1/R ~ fj, c± 15 MeV and M s = 1 TeV this 
square is given by 



1(T 2 



N*N ~ 1 + 0.05 (54) 



It constraints r to be in a subset of the range (||o[) which, however, is still comfortably 
large. 

Now that we have basically fixed all the lepton number conserving quantities in our 
model we should analyze what we need to do to generate small masses for the neutrinos. In 
the previous section we have shown that we can generate neutrino masses by introducing a 
lepton number breaking bulk Majorana mass M and lepton number violating brane-bulk 
masses m±. Both mass terms violate lepton number by two units. For the simple case 
of one neutrino flavor the resulting neutrino mass has been given in eq. (p9|). Using our 



favorite values 1/R ~ \x ~ 15 MeV and rh± ~ rh T where rh T has been specified in eq. (51) 
we find 

/10~ 2 \ f rh T M\ , 
m ~ f -y- J f -2-r 1 + — J (1.5-2.0) MeV . (55) 

Note that the two ratios rh T /rh T and M/fj, of lepton number violating and conserving 
quantities on the brane and in the bulk, respectively, enter with the same strength. This fits 
nicely to the fact that they both break lepton number by two units. More quantitatively, 
to have a neutrino mass ttiq < 1 eV we need both ratios to be less than roughly 10 _4 r. 
Such a number may, for example, arise from a suppression of the form < (f> > /M 8 where 
(p is a boson that carries lepton number and takes a VEV of the order, say 1/R. 

Although the main focus of the paper is on the particle physics aspects of our model, 
we would like to briefly address constraints from big-bang nucleosynthesis and astro- 
physics. In the early universe, the Kaluza-Klein mode with mode number n and mixing 
Ui n decouples at a temperature of 

T deCjn ~ lMeV|f/ m |- 2 / 3 . (56) 

From eq. (|51~1 ) the largest mixing angles are those for i = r. Using these angles and 
fj, = 1/R ~ 15 MeV as above we find 

1 /3 

T d ec,n ^ (10-11) MeV (^) (1 + n 2 ) 1 / 3 . (57) 

It follows than the decoupling temperature is lower than the mass for all Kaluza-Klein 
modes. In particular, the decoupling temperature increases with a smaller power of the 
mode number n than the masses M n ~ 15 M.eVy/1 + n 2 of the Kaluza-Klein modes. 
As a consequence, modes with large n decouple when they are highly non-relativistic 
and are strongly diluted. Unfortunately, this suppression is not strong enough for low 
n modes. Therefore, we have to assume that the bulk is empty at some temperature 



T < 7d CC] o — 10 MeV as it is customary for models with large additional dimensions [21]. 
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Taking this temperature sufficiently low (but above 1 MeV, of course) recreation of Kaluza- 
Klein particles is Boltzmann-suppressed. As a consequence, their thermal production rate 
can always be kept below the Hubble rate. 

Supernova cooling provides another constraint on our model. For the single-particle 
explanation of the KARMEN anomaly the supernova energy loss induced by the X particle 
has been estimated in ref. [p2j] . From this estimate, it is concluded there, that the single 
particle explanation of the KARMEN anomaly is ruled out. Our model is even more 
problematic since we have a tower of particles rather than just a single one. However, a 



closer analysis of the supernova energy loss [23| is likely to weaken the bound given in 
ref. p^j. This together with the assumption of a moderately lower supernova temperature 
might render our model consistent with the supernova constraint. In any case, we believe 
that our model is of interest from the viewpoint of particle physics and should be compared 
with the direct experimental information on neutrino properties. 



6 Conclusions 

In this paper, we have shown that the KARMEN anomaly can be understood in the 
context of a six-dimensional brane-world model. The slow particle X, responsible for the 
anomaly, is identified with a specific Kaluza-Klein excitation of a bulk fermion that, from 
a four-dimensional point of view, appears as a sterile neutrino. We have pointed out that 
in any interpretation of the KARMEN anomaly based on a single slow particle produced 
in the ir + — * [i + X decay, the X mass is fine-tuned to the mass difference m n — 
with an accuracy of order 10 -4 . This means that 5mx/(m-K — Trip) ~ 1.8 10~ 4 , where 
5mx = (rn n — m^) — mx- Such a problem can be significantly alleviated, although only 
partially solved, in the approach we propose, since the X particle is just one of the many 
excitation of the bulk fermion. It is therefore relatively more probable that one out of the 
many states has a mass that falls into the critical range. 

The phenomenology of the model with respect to the KARMEN experiment depends 
on the average separation Atokk between two Kaluza-Klein states at the scale mx ■ In the 
limit where Awkk is large, that is AmxK = O (m % — mj, the model is indistinguishable 
from models with a single sterile neutrino. In this limit, although the fine-tuning is not 
improved, we still have a model that provides a theoretically well-motivated origin for 
the X particle. Moreover, as a quite non-trivial feature, the model allows to incorporate 
the large mass scale of the X particle as well as the small neutrino masses needed for the 
standard oscillation phenomenology. The limit of small AmKK, that is Awkk = O (Smx), 
is ruled out because it would lead to a modification of the KARMEN time spectrum in the 
region (0.6-3.1) [is. For intermediate values of Arn-KK; Ath-kk = O (1) MeV or more, the 
fine-tuning problem is significantly alleviated and the particles lighter than X could give 
a small contribution to the early KARMEN time spectrum. Therefore, a detailed study 
of the time distribution of events at early times after the end of the proton pulse could 
test the model if the mass separation AmKK happens not to be too large. An upgraded 
detector with tracking capability devoted to the study of the anomaly could certainly 
explore this possibility. Needless to say, the detection of a signal at early times would 
represent a strong hint for the existence of a brane-world. 

On the model building side, we have identified two possible structures for the sizes 
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of the additional dimensions. The first case represents an effectively five-dimensional 
model with intermediate string scale and the remaining five dimensions being close to 
the string scale. The second one has a TeV string scale, two large dimensions in the 
MeV -1 range and the scale of the remaining four dimensions being two to three orders of 
magnitude larger. Therefore, in this case the space-time is effectively six-dimensional in an 
intermediate energy range. We have described in detail a six-dimensional example. As in 
all brane-world neutrino models, we have considered bulk neutrinos whose Lorenz invariant 
Dirac and Majorana masses are suppressed relative to the string scale. In particular, we 
have taken the magnitude of the Dirac mass term to be of the order of the inverse size 
of the additional dimensions and we have further suppressed the Majorana mass term by 
means of a U(l) symmetry (lepton number). Despite the introduction of the bulk Dirac 
mass term, we have found that the model accommodates three light Majorana neutrinos 
which would be massless in the limit of unbroken U(l) symmetry. The light states are 
predominantly given by the electroweak eigenstates with small but sizeable admixture of 
Kaluza-Klein modes. We have therefore attributed the solar and atmospheric neutrino 
phenomenology to those light states. On the other hand the small heavy component of 
the flavor eigenstates accounts for the KARMEN anomaly. Therefore, the introduction 
of the Dirac mass term makes our model rather different from the models considered in 
the literature so far, where the massless modes either do not exist or decouple from the 
left-handed neutrinos. Moreover, the neutrino masses and the mixing of neutrinos and 
bulk states are controlled by different parameters. We believe that such a scheme can be 
of general interest for brane-world models of neutrino masses. 
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